New exact solutions of the standard pairing model for well-deformed nuclei 
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A new step-by-step diagonalization procedure for evaluating exact solutions of the nuclear de- 
formed mean-field plus pairing interaction model is proposed via a simple Bethe ansatz in each step 
from which the eigenvalues and corresponding eigenstates can be obtained progressively. This new 
approach draws upon an observation that the original one- plus two-body problem in a fc-particle 
Hilbert subspace can be mapped unto a one-body grand hard-core boson picture that can be solved 
step by step with a simple Bethe ansatz known from earlier work. Based on this new procedure, it 
is further shown that the extended pairing model for deformed nuclei [Phys. Rev. Lett. 92, 112503 
(2004) ] is similar to the standard pairing model with the first step approximation, in which only 
the lowest energy eigenstate of the standard pure pairing interaction part is taken into considera- 
tion. Our analysis show that the standard pairing model with the first step approximation displays 
similar pair structures of first few exact low-lying states of the model, which, therefore, provides a 
link between the two models. 
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I. INTRODUCTION 

Pairing is an important residual interaction when a mean-field approach is used as a starting approximation for a 
description of nuclear structure. In particular, pairing correlations are essential for a description of binding energies, 
odd-even effects, single-particle occupancies, excitation spectra, electromagnetic transition rates, beta-decay proba- 
^ ' bilities, transfer reaction amplitudes, low-lying collective modes, level densities, and moments of inertia, etcJ 1_3 J Two 
commonly used methods, the Bardeen-Cooper-Schrieffer (BCS) [H and Hartree-Fock-Bogolyubov (HFB) 0] tech- 
nique for finding approximate solutions are well known. The limitations of these methods, when applied in nuclear 
] physics, are well understood, I 5 1 which is also the case when using these methods to determine the energy spectra 
J> ■ of nano-scale metallic grains. I 6-7 ' Various procedures have been used to correct the approximation deficiencies, such 
as particle-number projected mean-field treatments, I 8-10 ' the use of coherent states, I 11 ' stochastic number projection 
techniques, I 12 ' statistical descriptions, I 13 ' treatments of the residual parts of the Hamiltonian in the random phase 
approximation, [ 14 ~ 15 1 and various recursive approaches. I 16-17 ' Typically, however, these procedures have been found 
only limited applicability because the results usually yield insufficient accuracy. Other methods are likewise limited 
by their own sets of complications. 

On the other hand, an exact treatment of the nuclear mean-field plus pairing type Hamiltonian was initiated 
by Richardson and Gaudin, known as the Richardson-Gaudin method J 18-21 ' Recently, extensions to the Richardson- 
' Gaudin theory have been made using the Bethe ansatz methodology, I 22-31 ' especially, an application of the Richardson- 
. Gaudin solution to Sm isotopes with less than seven pairs of valence nucleons was made. I 32 ' Though these approaches 
show that the mean-field plus pairing model is exactly solvable, the solutions are generally not simple and normally 
require extensive numerical work, especially when the number of levels and valence pairs are large in spite of the 
recent efforts in improving the procedure J 33 ' In [34], an extended pairing model with many-pair interaction terms was 
proposed, which can be solved based on a simpler Bethe ansatz, and describes even-odd mass differences in 154 - 171 Yb 
isotopes rather well. However, it was not clear why the paring interaction strength in the model should be drastically 
reduced with increasing of the number of valence nucleon pairs in order to fit experimental data of even-odd mass 
differences. Moreover, since the standard form of pairing Hamiltonian is commonly adopted, it should be interesting 
to see whether the extended paring model in some aspects covers main features of the standard pairing model. 

In the following, a new step-by-step diagonalization procedure for evaluating exact solutions of the nuclear 
deformed mean-field plus pairing interaction model will be proposed via a simple Bethe ansatz in each step from 
which the eigenvalues and corresponding eigenstates can be obtained progressively, which is shown in Sec. II. In Sec. 
Ill, as a demonstration of the procedure, a system with p = 10 levels will be analyzed. In Sec. IV, it will be shown 
that the extended pairing model for deformed nuclei is similar to the standard pairing model with the first step 
approximation, in which only the lowest energy eigenstate of the standard pure pairing interaction part is taken into 
consideration. Our analysis reveals that the standard pairing Hamiltonian with the first step approximation displays 
similar pair structures of first few exact low-lying states of the model, which, therefore, provides a link between the 
two models. A short discussion regarding implications of our results is given in Sec. V. 



X 

5-H 



2 



II. A NEW DIAGONALIZATION PROCEDURE FOR THE MEAN-FIELD PLUS STANDARD 

PAIRING HAMILTONIAN 

The deformed mean-field plus standard pairing model Hamiltonian is given by 

H = E e A - GS+S- = ]T ejnj - G £ , (1) 

where p is the total number of levels considered, G > is the overall pairing strength, {e.,} are single-particle 
energies taken from any deformed mean-field, such as the Nilsson model! 2 ! or th e relativistic mean-field theory, I 35-37 ! 
hj = Cj^Cji + c j| c il is the fermion number operator for the j-th level, and = cj^ct^ (S^ = (S^Y = c^Cii) are pair 
creation (annihilation) operators. The up and down arrows in these expressions refer to time-reversed states. Since 
the formalism for even-odd systems is similar, in the following, we only focus on the even-even seniority zero case. 

Since each deformed level can be occupied by no more than a single pair due to the Pauli principle, Hamiltonian 
(1) is also equivalent to a finite-site hard-core Bose-Hubbard model with infinite range hopping and infinite on-site 
repulsion. Let Bf ii2 ... ik — 5*+ 5+ ■ ■ ■ with 1 < i\ < i 2 < ■ ■ ■ < ik < V- I n the following, we set p — (i\i 2 ■ ■■ik) to 
be the /U-th normal order sequence with 1 < i\ < i 2 < ■ ■ ■ < ik < P- The operator can be regarded as a grand 
hard-core boson creation operator J 34 l The total number of these operators is p\/((p — k)\k\) which exactly equals the 
dimension of the Hilbert subspace of k pairs with no double occupancy. 

For fc-pair excitations, by using the standard second quantization formalism, the Hamiltonian (1) can effectively be 
reduced to the following 'one-body' Hamiltonian in the grand hard-core boson picture: 

H k =J2^\i^»B+B v - GJ2^\S + S-W)B+B^ (2) 

liv j=l p,u 

where B^ = (-B+)t, because any fc-pair eigenstate of (1) can be expanded in terms of single grand hard-core boson 
states = £> + |0)}, where |0) is the pair vacuum state. In (2), the mean-field one-body term in the Hilbert subspace 
spanned by {|^) = B+\0)} is diagonal with the matrix elements 

p k 

while the matrix elements of the pairing interaction term are 

= £>|p/2 - q P k) QQ (p/2 - q pk\S+S-\p/2 - q pk) Q Q {p/2 - q pk\v), (4) 
<?p 

in which 

Q (p/2 - q pk\S+S-\p/2 - q pk) Q =h[ q) =(k-q)(p-k-q+ 1) (5) 

is matrix element of S + S~ in the Racah quasi-spin formalism! 38 ] with the total quasi-spin 5q = p/2 — q, where 
q = 0, 1, • • • , min[fc,p — k], and p is an additional quantum number needed in distinguishing from different states with 
the same quasi-spin for given p and k. For given p, the total number of different quasi-spin states with the same 
Sq = p/2 — q is given byl 39 l 

(p - 2q + l)p\ 

^ ( P -q+l){p-q)\qV U 

Furthermore, aff — (p\p/2 — q pk)Q used in (4) is an overlap of the quasi-spin state \p/2 — q pk)Q with p-th single 
grand hard-core boson state \p), which can be chosen as real. It can easily be verified that the total number of different 
quasi-spin states equals exactly to the dimension of the Hilbert subspace of k pairs with no double occupancy: 
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min[k,p— k] . 



Thus, (2) can explicitly be written as 



d min[fc,p— k] uj q 

H k = ^2e,B+B,-G £ ^g^a^B+B,. (8) 

/j=1 g=0 p=l A 1 ^ 

In order to diagonalizc the Hamiltonian (8), let us consider a simpler Hamiltonian with only first term and q = 
part in the second term of (8): 

d 

h a = J2 XpB+Bp - Ghf afafB+B u . (9) 
As shown in [34, 40], The Hamiltonian (9) can be digonalized into the following form: 

h =J2 E {t ^D + {E {tq ^)D{E {t ^) (10) 

TO 

with 



D+(E™) = W-L Y ^-^B+ (11) 



where E^ is the ro-th root of the following equation: 



and 7V T0 is the normalization constant obtained from 



(12) 



K 1 ) 2 



= Kt>K - (13) 



In this case, (12) will provide with exactly d different roots E^ as long as all combinations of the single-particle 
energies e M = X)t=i e n are different for all fc-pair excitation cases. Fortunately, this is always the case when single- 
particle energies {e^} are generated from any deformed mean-field theory. 

Since (8) and (9) should be diagonalized within the single grand hardcore boson subspace spanned by {-B+|0)}, the 
effective commutation relations needed to prove that (9) can indeed be expressed in the form shown in (10) are 

[B v ,B+]=5 liV , (14) 
which are only valid when they are applied onto the vacuum state. Using (9) and the ansatz (11), we have 



[h ,D+(E^)] = E^ D +(E^) Ghf £ J^ {to) ) E <W (15) 
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Though a direct proof is in demand, it can be checked numerically with any set of parameters, Gh® \ {a^ 1 }, and 
{e^} with ei 7^ e 2 • • • =/= Id, that the orthnormal condition (13) is automatically satisfied when i?( T °) satisfies (12). 
Therefore, (9) can indeed be expressed as that shown in (10) as long as (12) is satisfied. 

In order to simplify our expression, in the following, the indices (q,p) are relabeled by r with r = (q,p). Thus, the 
Hamiltonian (8) can be rewritten as 

d-l 

4 = ^£Wfl+(£W )J)( £( T .)) _ G ^4^ a ;a^;B, (16) 

to r=l ilv 

Then, we similarly have 

J2e^D+{E^)D{E^) - Gh^^alalB+Bv = 

TO fJ-V 

J2e (to) d + (e^)d(e^) -Gh { k 1] J2 &t a t >d + (e^)d{e^) = J2 e(ti)d+ ( e(ti) ) d ( e(ti) )> 

T T T Q T X 

where 

Ato = V^^ 2 ^-^ o) ' (18) 

V l-i 



V T 

E( Tl *> is the n-th root of the following equation: 

GhPY , {Ato)2 ( , =1, (20) 

TO 

and Af T1 is the normalization constant obtained from 



V ( At °) 2 — a ,\r foil 

^ (E^ - E^))(E^') - £«)) "-nr^n- 
Thus, using the results shown in (17)-(21) and following the above procedure consecutively, we finally have 

d-l 

H k = J2 E(ro)D+ ( E{To) ) D ( EiTo) ) - G Y, h k } Y, ar ^ B ^ Bv = E^-^D+(E { - Td -^)D(E { - Td -^), (22) 

where 



Td — 2 



with 



, _ /j_ v ff A - v f 24 ) 

' - - \r_ 2^ 11 R (t„) _ r(t„ +1 ) 2^ of. - rto ' ^ Z4 J 



Ts V A/V ^ 11 S< T -) - E'K+i) ^ 2e„ - ' 

s ToTl---T s _l !/ = ^ 
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c+(£ """ ) ^vS? ^'-^" , I ' t(E ' r '' ) ' (25) 



and 



^ (£('.) - £(T.+i)) T [i(T.) _ B (ri +1 )) 5 ^+iri +1 ^r. +1 , (27) 

for s = 0, 1, 2, • ■ • , d— 2. Hence, after d steps, the Hamiltonian (8) is diagonalized as shown in (22), of which eigenstate 



is 



|fc,r d _ 1 )=^+(£;^- 1 ))|0) (28) 

with the corresponding eigen-energy E^ Ti -^ . 

This new step-by-step diagonalization procedure needs at most d steps to get final exact results, but in each 
step the corresponding Bcthe ansatz equation (26) contains only one variable, of which roots can easily be obtained 
numerically similar to what is required in the extended pairing model proposed previously, I 34 ' and in the TDA and 
RPA approximations with separable potentials.! 2 ] Though this method may be unpractical for large size systems 
because one needs to get all d roots from the equation (26) in each step, this procedure can also be used to check 
contributions from pairing potential in the Racah quasi-spin formalism for different q of the second term in (8), and 
is certainly applicable to relatively small systems. Actually, for k pair excitation, though each term with different q 
from the second term of (8) will contribute to the final eigen-energy and correlate with eigenstates, the first few of 
these terms are key to determine properties of the first few low-lying states of the model as will be shown in the next 
section. 

III. A NUMERICAL EXAMPLE FOR p = 10 

In this section, we will apply this new step-by-step diagonalization procedure to the deformed mean-field plus 
standard pairing model for p = 10 levels with number of pairs k = 1, 2, • • • , 10, in which the single particle energies 
are given by ej = i + Xi for i = 1, 2, ■ • • , 10, where \i are random numbers within the interval (0, 1) to avoid accidental 
degeneracy, and the pairing strength is set to be G — 0.5. since > > ■ ■ ■ > h d ~ x is always satisfied, the lowest 
quasi-spin term with q = from the pairing potential should be most important to the first few of eigenstates of the 
model, which is indeed the case as can be seen from results shown in Table I. 

In each step, we need to calculate the overlaps of the quasi-spin states \p/2 — q pk)q with /x-th single grand hard-core 
boson states \p). For q = 0, the fc-pair state 



\ P /2 k) Q = y/(p ky./(k\p\)(s + ) k \o) = y/k\(p - k)\/ P \ Y, s+sf-.-sl\o) = ^T/dJ2B+\o) (29) 

l<il<i-2<---<ik<P A 1 

is the eigen-state of the operator S + S~ . Thus we have 



af = (»\p/2k) Q = ^TJd. (30) 

For q > 1, the quasi-spin states \p/2 — q pk)Q can be obtained by directly diagonalizing S + S~ as shown in (5), or 
by using representation theory of SUq(2) x S p summarized in [39]. Then, one can use them to calculate the overlaps 
oflP = (Mb/2 - q pk) Q . 

In Table I, we list first 5 eigenenergies and overlaps of the eigenstates with the corresponding exact ones for number 
of pairs k — 1, 2, • • • , 10 calculated with only h term involved, which is called the first step approximation. With 
the first step approximation, it can be seen from Table I that results of k — 1 and k — 10 cases are exact because 
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TABLE I: First 5 eigenenergies of the standard pairing model with p = 10 levels for k = 0, 1, 2 • • • , 10 obtained from the first 
step approximation (appro.) and compared with the corresponding exact results (exact), and overlap (olp) of the eigenstates 
obtained from the first step approximation with the corresponding exact ones, in which the single-particle energies ei = 1.706, 
e 2 = 2.754, e 3 = 3.440, e 4 = 4.349, e 5 = 5.743, e 6 = 6.604, e 7 = 7.591, e 8 = 8.959, e 9 = 9.335, e 10 = 10.125, and pairing strength 
G — 0.5, where the single-particle energies and G are given in arbitrary units. 



k=l 






k=2 




k=3 




k=4 






k=5 




exact appro. 


olp 


exact 


appro. 


olp 


exact appro, olp 


exact 


appro. 


olp 


exact 


appro. 


olp 


Ei 2.255 2.255 


100% 


6.662 


8.133 


98% 


12.873 15.423 93% 


21.101 


24.272 


91% 


31.856 


35.748 


88% 


E 2 4.797 4.797 


100% 


8.999 


9.808 


97% 


15.508 17.213 92% 


24.293 


26.909 


90% 


34.529 


37.447 


88% 


E z 6.438 6.438 


100% 


10.756 


11.392 


78% 


16.961 18.576 89% 


26.161 


28.312 


65% 


36.407 


39.209 


83% 


Ei 8.272 8.272 


100% 


13.128 


12.2618 75% 


17.973 19.832 82% 


26.176 


29.058 


67% 


36.961 


40.167 


81% 


E 5 10.969 10.969 100% 


13.375 


13.625 


85% 


18.936 20.598 65% 


27.621 


29.882 


62% 


38.828 


41.206 


55% 


k=6 






k=7 




k=8 




k=9 






k=10 




exact appro. 


olp 


exact 


appro. 


olp 


exact appro, olp 


exact 


appro. 


olp 


exact 


appro. 


olp 


Si 44.638 48.868 


91% 


59.532 


63.713 


96% 


76.780 79.982 96% 


95.625 


97.091 


95% 


116.212 116.212 100% 


E 2 47.415 50.152 


90% 


62.610 


65.911 


93% 


78.949 82.641 97% 


97.888 


101.652 96% 








E 3 49.165 52.261 


81% 


63.852 


67.440 


82% 


80.277 83.861 81% 


99.015 


102.961 99% 








Ei 50.908 53.366 


76% 


64.573 


69.232 


71% 


81.182 85.191 82% 


101.540 105.221 98% 








E 5 51.728 54.227 


74% 


65.263 


69.651 


55% 


83.148 86.404 75% 


103.552 


107.341 98% 









TABLE II: First 5 eigenenergies of the standard pairing model with p = 10 levels for k = 5 obtained from step-by-step 
diagonalization procedure, in which the parameters used are the same as those shown in Table I, where only with q = 
0, 1, 2, • • • , s terms from the pairing potential are involved in the s-th step approximation as described in the previous section. 





1st step 




2nd step 




3rd step 




4th step 




exact 






q=0 




q=0,l 




q=0,l,2 




q=0,l,2,3 




q=0,l,2,3,4 






eigenvalue 


olp 


eigenvalue 


olp 


eigenvalue 


olp 


eigenvalue 


olp 


eigenvalue 


olp 




35.748 


88% 


33.862 


98% 


32.481 


99% 


32.176 


99% 


31.856 


100% 


Ei 


37.447 


88% 


36.566 


97% 


35.631 


99% 


35.102 


99% 


34.528 


100% 


Ez 


39.209 


83% 


38.398 


90% 


37.451 


99% 


36.847 


99% 


36.407 


100% 


Ei 


40.167 


81% 


39.422 


89% 


37.942 


99% 


37.637 


99% 


36.961 


100% 


E 5 


41.206 


55% 


40.821 


57% 


39.779 


97% 


38.819 


97% 


38.828 


100% 



/j(<j) — for q > 1. The approximate energy eigenvalues will gradually greater than the corresponding exact ones with 
increasing of the number of pairs k since pairing potential terms with q = 1, 2, • ■ • , min[p — k, k] will contribute 
more and more to the final eigenenergies. However, the overlaps of the eigenstates with the corresponding exact ones 
are always greater than 88% for the ground and the first excited states for any number of pairs k. Therefore, 
term of the pairing potential is dominant in determining pairing structure of the first two excitation states in the 
model though the corresponding energy eigenvalues arc different from the exact ones. 

Since the largest deviation of the energy eigenvalues from the exact ones occurs at the half-filling case, using 
the procedure shown in the previous section, we calculated energy eigenvalues step by step for k — 5 case with 
q = 0, 1, • • • , 4 since — 0, of which the results are shown in Tabic II. For given q, there are actually uj q sub-steps 
involved in the diagonalization process according to the procedure shown in the previous section. It can be seen 
from Tabic II that the overlaps of the first 4 eigenstates with the corresponding exact ones will reach 99% after three 
diagonalization steps though there is still deviation in eigenenergies, which shows that the first few terms with 
q = 0, 1, 2 are key to determine pair structure of the first few low-lying states of the model. While high- lying quasi-spin 
states mainly correlate with high excited states of the model and keep low part of the spectrum less affected. 

IV. COMPARISON WITH THE EXTENDED PAIRING MODEL 

For k pair excitations, if only term from the standard pairing potential is considered for the standard pairing 
Hamiltonian (8), namely 
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= £ ^B+B, - G(k(p - k + l)/d) £ B+B v , (31) 



following (9)-(13), eigenstate of (31) cen be written as 

d 



^) = E^r^F)^ +|0) ' (32) 



where r is an additional quantum number used to distinguish different excitation states, and E^ is an unknown 
variable to be determined in diagonalizing (31). In solving the following eigen-equation 

H^\k;r)=E^\k;r) (33) 
with BpB+ 10) = <5 M „|0), it shows that the variable E^ must satisfy the following equation: 

d i 

_ =1, (34) 

and the eigen-energy E^ = E^ T \ Thus, the additional quantum number r labels different roots of (34). This is 
the so-called the first step approximation shown in section II. The solution is complete so long as all combinations of 
the single-particle energies J2t=i e n arc different for all fc-pair excitation cases. Fortunately, this is always the case 
when single-particle energies {ej} are generated from any deformed mean- field theory. Since the single grand particle 
energies 2e^ are all different, there are exactly p\/((p — k)\k\) distinct roots in (34). The resultant eigenstates (32), 
which are mutually orthogonal but not normalized, satisfy 

(fc; 7~|fc; t') = 5 TT 'Af T , (35) 

where 



It follows that the normalized eigenstate can be expressed as |fc;r) = \Jl/N T \k; t). 
As shown in [34], a Nilsson mean-field plus extended pairing interaction Hamiltonian 



H c . = £ e jnj - G cx £ St Si G cx ]T £ S £SZ ■ ■ ■ St p S lp+l S lp+2 ■ ■ ■ a l2p , (37) 

where no pair of indices among {ii, ■ ■ ■ , i2p\ are the same for any p, can also be solved exactly by using a simple 
Bethe ansatz that is similar to what is proposed in this work. Besides the usual Nilsson mean-field and the standard 
pairing interaction, this form includes many-pair hopping terms that allow nucleon pairs to simultaneously scatter 
(hop) between and among different Nilsson levels. Furthermore, the extended pairing interaction Hamiltonian (37) 
can be used to describe even-odd mass differences rather well as long as the extended pairing interaction strength 
G e x decreases with an increasing number of pairs k. It follows from this that it is interesting to compare results of 
the deformed mean-field plus standard pairing Hamiltonian (1) with those from the extended pairing modelJ 34 ' And 
indeed, it is not difficult to show that the expressions for eigenstates of the extended pairing Hamiltonian and those 
of the standard pairing Hamiltonian (1) in the first step approximation are the same. For fc-pair excitations, the 
eigenenergies E^\ex) of the extended pairing Hamiltonian (37) arc given by 

4 T) (ex)= E M-(k-l)G cx , (38) 
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where i?J x is the r-th root of the Bethe ansatz equation, 

Ge*£— -^y = l. (39) 

A comparison of (39) with the Bcthc ansatz equation (32) for the standard pairing Hamiltonian in the first step 
approximation (31) shows that the two Hamiltonians yield exactly the same excitation energies and the corresponding 
eigenstates so long as the parameter G ex in the extended pairing model and the parameter G in the standard pairing 
Hamiltonian (1) satisfy the following relation: 

G cx = {{ P -ky.k\(p-k + l)k/p\)G. (40) 

Furthermore, while the ground states of the two Hamiltonians are also the same, the ground-state energies arc different. 
However, once the overall pairing strength G is fixed, and the parameter G ex is chosen according to (40), it is easy 
to show that the difference between the ground-state energy of the extended model and that of the standard pairing 
model in the first step approximation is given by 

4 g) (ex)-4 s > = -(fc-l)G 0X . (41) 

This expression shows that the extended pairing interaction contributes a little more attraction among valence pairs 
than the standard pairing interaction in the first step approximation, but reproduces excitation energies exactly the 
same as those in the standard pairing model with first step approximation. Since G cx decreases drastically with 
increasing of k toward the half-filling, the ground sate energy difference of the two Hamiltonians becomes negligible 
with increasing number of pairs k with k < [p/2] — 1 when p is even, and k < [p/2] when p is odd, where [x] denotes 
integer part of x. 

As an example, the ground-state energy difference (41) of the two Hamiltonians in the sixth (82-126) major shell 
with the standard pairing strength G = 0.2MeV, which is a typical parameter value for describing deformed nuclei in 
this region, shows that the ground-state energy difference of the two Hamiltonians are rather small in this case. The 
largest deviation of the ground-state energy of the two Hamiltonians is at k = 2 with s[ g ^(ex) — = — 36.3636keV. 
Notwithstanding, since the only difference between the two Hamiltonians, so long as G ox is taken to be related to G 
by prescription (40) , is in the overall binding energy, and since an analytic expression for this difference in also known 
in terms of G ex through (41), for practical purposes the two Hamiltonians yield the same results, even though the 
Hamiltonians are quite different. This in itself is interesting, since it shows that a many-pair interaction Hamiltonian 
can have identical solutions to the two-pair interaction with truncations. Obviously, it follows that for such systems 
the structure of fixed-Z (isotopic) and fixed-N (isotonic) chains follow solely from the structure the simplest single-pair 
member of the chain and simple "pair-counting" factors related to the pairing interaction strength and single-particle 
energies. 

Thus, we conclude that, basically, the extended pairing model is different from the standard pairing model. However, 
if only a first few eigenstates are considered, the pair structure of these states in the two models are similar, especially 
in ground state, as can be seen from analysis of the overlaps in the previous section. It can be expected that the 
difference of the two models will be negligible when number of pairs k or pairing interaction strength G is small. In 
addition, since the extended pairing model can be solved exactly with a single one variable equation (39), which is 
simpler than the Richardson- Gaudin equations with k variables for the standard pairing Hamiltonian, the extended 
pairing model can be applied to relative large systems, especially when one only want to know a first few low-lying 
eigenstates and corresponding eigenenergies. 

V. CONCLUSION 

A new step-by-step diagonalization procedure for evaluating exact solutions of the nuclear deformed mean-field 
plus pairing interaction model is proposed via a simple Bethe ansatz in each step from which the eigenvalues and 
corresponding eigenstates can be obtained progressively. This new approach draws upon an observation that the 
original one- plus two-body problem in a fc-particlc Hilbcrt subspace can be mapped unto a one-body grand hard-core 
boson picture that can be solved step by step with a simple Bethe ansatz known from earlier work, in which one only 
needs to solve a single variable non-linear equation instead of a set of coupled non-linear equations with k variables as 
is required, for example, within the framework of the well-known Richardson-Gaudin method. Though this method 
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may be unpractical for large size systems because one needs to get all d roots from the Bethe ansatz equation in each 
step, this procedure can be used to check contributions from pairing potential in the Racah quasi-spin formalism, and 
is certainly applicable to relatively small systems. 

As is shown in the example with p = 10 levels, though each term with different q from the pure pairing interaction 
will contribute to the final eigen-energy and correlate with eigenstates, a first few of these terms are key to determine 
a first few low-lying states of the model. While high-lying quasi-spin states mainly correlate with high excited states 
of the model and keep low part of the spectrum less affected. 

Based on this new procedure, it is further shown that the extended pairing model for deformed nuclei! 34 ' is similar 
to the standard pairing model with the first step approximation, in which only the lowest energy eigenstate of the 
standard pure pairing interaction part is taken into consideration. Our analysis show that the standard pairing 
Hamiltonian with the first step approximation displays similar pair structures of a first few low-lying states of the 
standard pairing model, which, therefore, provides a link between the two models. 

Furthermore, the new method proposed is not limited to the deformed mean-field plus pairing problem only, as it 
should also prove useful for solving a much large class of quantum many-body problems in which model Hamiltonians 
are described by 

H = H + XH U (42) 

where A is a real parameter, and H and Hi do not commute, [Ho, Hi] ^ 0. According to our procedure, if the 
particle number is a conserved quantity, and Ho and Hi can be diagonalized independently in a fc-particlc basis, 
then the Hamiltonian (42) is exactly solvable by using the step-by-step exact diagonalization procedure. Moreover, 
the method can also be extended to deal with Hamiltonians with more than two non-commutative terms by using a 
similar procedure consecutively. Research in this direction is in progress. 
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